Non-Abelian gauge fields and topological insulators in shaken optical lattices 
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Time-periodic driving like lattice shaking offers a low-demanding method to generate artificial 
gauge fields in optical lattices. We identify the relevant symmetries that have to be broken by the 
driving function for that purpose and demonstrate the power of this method by making concrete 
proposals for its application to two-dimensional lattice systems: We show how to tune frustration 
and how to create and control band touching points like Dirac cones in the shaken kagome lattice. 
We propose the realization of a topological and a quantum spin Hall insulator in a shaken spin- 
dependent hexagonal lattice. We describe how strong artificial magnetic fields can be achieved for 
example in a square lattice by employing superlattice modulation. Finally, exemplified on a shaken 
spin-dependent square lattice, we develop a method to create strong non-Abelian gauge fields. 



Topological order and topological insulators [lj are cur- 
rently in the center of interest of quantum physics, espe- 
cially because of their possible applications in quantum 
information and spintronics [2 J. For this reason, there 
is an ongoing search for feasible realizations of such sys- 
tems in- and outside of solid-state physics. Here, ultra- 
cold ground-state atoms provide a promising playground 
[3j (although Rydberg-excited atoms [4], trapped ions 
[5J, and photons in nano-structured materials [6j offer 
interesting alternatives). Typically, topological effects 
require ultra-strong gauge fields or spin-orbit-like cou- 
plings. There are several ways to achieve these with 
ultracold atoms, from trap rotation [7J, microrotation 
[8], to Berry phase imprinting |9J. In optical lattices, 
combining laser-induced tunneling with superlattice tech- 
niques allows for strong Abelian [10J and non-Abelian 
[TT] gauge-fields and for the realization of topological in- 
sulators [12]. So far, the first lattice experiments led to 
the creation of staggered flux lattices [13J. Many other 
groups follow this direction of research [14J. 

Recently, there has been a burst of interest in another, 
experimentally less demanding, approach, namely peri- 
odic lattice shaking. Sinusoidal shaking leads to a change 
of strength, or even sign of the tunneling and allows to 
control the Mott-insulator-supernuid transition [15] [16] 
(for a recent work in hexagonal geometry, see [17J). While 
in the square lattice this introduces neither frustration 
nor synthetic gauge fields, in the triangular lattice a sign- 
change of the tunneling is equivalent to a 7r-flux Abelian 
field [18J. Such a system mimics frustrated antiferro- 
magnetism, classical for weakly interacting bosons [19J, 
and quantum in the hard-core boson limit [20J, where 
it is expected to exhibit exotic spin-liquid phases [2T] . 
Recently, it was demonstrated that by breaking tempo- 
ral symmetries of the shaking trajectory, one can cre- 
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ate phases of the tunneling in an optical lattice [221 [23] 
(see also Ref. [24J), and that in this way tunable Abelian 
fluxes through triangular plaquettes may be generated 
|22j . In this letter, we discuss non-trivial generalizations 
of this approach that involve also AC-induced tunneling 
and spinful particles. This allows us to simulate Abelian 
and even non-Abelian SU(2) gauge-fields in different lat- 
tice geometries, as well as topological insulators. To this, 
we employ non-standard optical lattices, like kagome and 
spin-dependent square and hexagonal lattices, and con- 
sider scenarios based on superlattice modulation. 

Basic scheme, and temporal symmetries. We con- 
sider a system of ultracold atoms in a driven optical 
lattice described by the Hubbard Hamiltonian H(t) = 

-J2(ij) Jij^l^j + Hi v i( t )^i + ^os with (bare) tunnel- 
ing matrix elements and annihilation and number 
operators a i and hi for particles (bosons or fermions) 
at site i\ H os collects on-site terms describing interac- 
tions or a weak static potential. The potential Vi(t) = 
vf(t) + VihhJ consists of two parts: a rapid periodic drive 
v?(t) = vf(t + T) of frequency u = 2n/T and zero 

time average (vf(t))T = with (-)t = \ Jq -dt; and 
(unlike in Ref. [22]) strong static energy offsets Viftjuo 
with integers V{. For huj ^> a large energy differ- 
ence Vijfiuj ^ (here and below we use the double-index 
shorthand X{ — Xj = Xij) practically prohibits tunnel- 
ing between i and j, unless the resonant periodic driving 
leads to AC-induced tunneling (ACT) [25], as it has been 
observed in recent experiments [26J. Later on, we will 
augment our model by a spin degree of freedom s =t, |. 

A gauge transformation U = exp (i^2 i Xi(*)^i)> where 
Xi(t) = Xi(t) ~ vM + 7 * with h X ?(t) = - f* drvf(r) + 
(Jq drvf(r))T and constants 7^, leads to the new Hamil- 
tonian H'(t) = U^HU— ihU^(dtU), which can be approx- 
imated by its time average H e ^ = — ^tf^l^j + ^os 
if tiuj is large compared to both the and the energy 
scales of H os . In this treatment, the initial energy offsets 
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Figure 1: (color online) (a-c) Lattice geometries involving 
triangular plaquettes pierced by an artificial magnetic flux 
4>v,a = ±0 (indicated by + and — ): (a) triangular lattice, 
(b) kagome lattice with tunneling J\ (J 2 ) along the darker 
horizontal (lighter diagonal) bonds, and (c) hexagonal lat- 
tice with nearest-neighbor ACT (solid lines) between shallow 
A- and deep B-sites and next-nearest-neighbor tunneling be- 
tween A-sites (dashed lines), (d) Driving function breaking 
symmetries (r) and (s). 

Viftu) enter via the effective tunneling matrix elements 
J if = Jij(e~^ Xij ^)T only, and in H e R all sites appear 
to have the same energy. In the undriven system, for 
Vij 7^ 0, the large energy difference Uijhuj suppresses tun- 
neling between sites i and j, and this fact is reflected in 
i^eff by a vanishing effective tunneling Jff = at van- 
ishing driving vfj =0. In turn, finite driving vfj ^ can 
establish coherent ACT with Jff ^ 0, where the energy 
difference v^hhj is bridged by quanta hu. 

The leitmotif of the present work is to use this control 
scheme to induce Peierls-type phases 

Oij = arg ((e-^x&W-^^+w))^ (1) 

that cannot be eliminated globally by choice of gauge, 
i.e., by adjusting the constants 7^. Such non-trivial 
phases correspond to artificial Abelian gauge fields; the 
gauge- invariant magnetic flux <pp £ (— 7r, 7r] piercing a 
lattice plaquette P is (modulo 2tt) obtained by summing 
the Oij around P. We find that the global reflection sym- 
metry (r) vf(—t — r) = vf(t — r) with respect to a global 
time r (using the choice 7^ = — z^cjt) implies trivial 
0^ = 0. Moreover, if ACT is not involved (z^j = 0), 
0^ = follows already from the local reflection sym- 
metry (r') vfj(—t — Tij) = vfj(t — Tij) with independent 
local times (since 7^ = VijUJT = 0, independent of r), 
or from the shift antisymmetry (s) vf(t— ^) = —vf(i) 
(choosing 7^ = 0) [27J. Therefore, ACT significantly re- 
duces the constraints on the driving function vf(t) for 
the creation of artificial magnetic fields. This is nicely 
exemplified by recent proposals where already simple si- 
nusoidal forcing [fulfilling (r') and (s)] leads to magnetic 
fields when combined with ACT - provided the temporal 
phase of the driving can be made site dependent [thus 
breaking (r)] [29J. In the following, we consider exper- 
imentally feasible scenarios where the whole system is 
driven in phase [such that both (r) and (s) are broken]. 

We will later generalize the scheme described in the 
two preceding paragraphs to the case of spin- 1/2 particles 
and show how non- Abelian gauge fields can be realized. 

Homogeneous forcing and triangular plaquettes. Let 
us consider a homogeneous time-periodic force F(t), such 



as an inertial force created by shaking the lattice along 
a periodic orbit. For Vi = 0, the driving potential 
vf(t) = —Vi'F(t) (with site position 7^) results in Peierls 
phases 0^ that only depend on the vector = Vi — re- 
connecting the two sites i and j, 0^ = f(rij). Using 
Eq. ([!]), one finds that /(— r^-) = —f(rij) and, therefore, 
homogeneous forcing cannot be used to create artificial 
magnetic fluxes through plaquettes with pairwise parallel 
edges. Since, however, generically 0^ depends in a non- 
linear fashion on [f(rij) is not of the form b-rij], one 
can use lattice shaking to induce a strong and tunable 
artificial magnetic flux 0v through, e.g., a downwards 
pointing triangular plaquette V. In the supplemental 
material [30J , we analytically compute this flux for unidi- 
rectional forcing. The inversion of the triangular plaque- 
tte V —> A reverses the sign of the flux, a = — </>v? such 
that staggered fluxes can be achieved in the triangular 
or kagome lattice as shown in Fig. [TJi and b. Since these 
flux configurations stem from homogeneous forcing they 
do not break the translational symmetry of the lattice. 

Tuning the staggered flux allows one to continuously 
control the degree of frustration in these lattices from 
none for zero-flux to maximum for 7r-flux [correspond- 
ing to ferromagnetic {—J if < 0) and antiferromagnetic 
coupling {—J if > 0), respectively]. The fully-frustrated 
regime gives rise to intriguing physics. For example, 
the flat lowest band of the kagome lattice makes the 
system extremely susceptible towards interaction-driven 
physics [31]; moreover, the case of hard-core bosons 
can be mapped to the spin-l/2-XY antiferromagnet [20J 
with possible spin-liquid ground states in the spatially 
anisotropic triangular lattice [21 J and still unexplored 
behavior in the kagome geometry. The ability to tune 
continuously between zero and maximum frustration de- 
scribed here can, thus, be a powerful tool for the adia- 
batic preparation of frustrated quantum phases. 

The realization of tunable staggered fluxes as shown in 
Fig. [T^i and b is also interesting in its own right. In the 
bosonic case, deviations from 7r-flux directly map to tun- 
able Dzyaloshinskii-Moriya couplings in the spin picture 
(see, e.g., [32J). Furthermore, for finite flux 0a = 0, the 
three bands of the kagome lattice feature a complex band- 
touching structure whose topology can be controlled by 
the driving. This is illustrated in Fig. |2^i for a lattice 
with I Jff\ equal to J\ {J 2) along the horizontal (other) 
bonds (see Fig. [TJd) . 

Topological and quantum spin Hall insulator. Such 
triangular plaquette fluxes can be used to engineer a 
topological insulator and a quantum spin Hall insulator. 
Consider a spin-dependent hexagonal optical lattice as 
sketched in Fig. where sites of the A (B) sublattice 
are energetically lifted (lowered) by AE/2 for t parti- 
cles, and vice versa for 4- particles [33J. Let us focus on 
non-interacting ^-particles first. For substantial detuning 
AE, we can assume that nearest-neighbor (NN) tunnel- 
ing (between A and B sites) is energetically suppressed 
and that next-NN (NNN) tunneling is relevant only be- 
tween sites of the "shallow" A sublattice. Now assume 





Figure 2: (color online) (a) The topology of band touching 
for the kagome lattice can be controlled by anisotropy J\ / J2 
and plaquette flux cj). The way and how often the three bands 
touch is depicted by the iconographic symbols, (b) Phase dia- 
gram of the hexagonal lattice as in Fig.[l]: with bare/undriven 
(next) nearest neighbor tunneling matrix elements J (J 7 ), 
subjected to a symmetry-breaking force of amplitude a and 
direction = cos((Pf)b x +sm((pF)B y , F(t). white: no Dirac 
points are present; gray: a small nearest-neighbor tunneling 
< 0.02 J renders the physics effectively ID. The colorbar en- 
codes the masses at two Dirac points, labeled as |m< | < |m> |. 
In the diagonally (horizontally) hatched region both masses 
are positive (negative). When the masses have opposite sign 
(un-hatched), the system is a topological insulator (or a quan- 
tum spin Hall insulator for two spin states). Inset: Position 
of Dirac points in k-space for LpF — 7r/30, indicating how they 
move and merge with a. 



that the system is driven resonantly by a time-periodic 
homogeneous force of frequency vab^ — (with inte- 
ger vab) that both establishes NN ACT and creates finite 
artificial fluxes through the triangular NNN plaquettes 
of the A sublattice ("+" and "-" in Fig. [TJ;). Introduc- 
ing Pauli matrices a for the sublattice degree of freedom, 
the effective Hamiltonian in momentum representation 
becomes # eff = Y.k°>kHk)a k where a ] k = (a^^k) 
and ft(fe) = 5ft (g(k)) a x -% (g(k)) <r y +g'(k) § (1 + s z <r z ). 
Here, s z = 1 and g^\h) = - ^(0 J$) exp(ifc-«W) with 
tfW denoting the three (six) vectors connecting an A site 
to its NN (NNN). Diagonalizing h(k) gives the dispersion 
relations e±(k) = \g'{k) ± yj\g{k)\ 2 + \g'(k)/2\ 2 for the 
two bands. 



Without NNN tunneling (g' = 0), the system can pos- 
sess a pair of band-touching points, i.e., g(ki^) = 0, with 
light-cone-like dispersion relation, so-called Dirac cones. 
A finite NNN g' (k) will split the bands at these points, 
and the Dirac- type dispersion relations found near 
acquire finite "masses" mi ,2 = g'{k\^)- If these have 
opposite sign, the lowest band possesses a finite Chern 
number (±1). Then, if the lowest band is entirely filled 
with t fermions, the system is a topological insulator with 
quantized Hall conductivity and robust chiral edge modes 
[34J (see also |35j). Repeating the above reasoning for \, 
particles, for which the role of A and B sites is inter- 
changed, one obtains the same result, but with s z = — 1 
and inverted Hall conductivity. Therefore, filling the low- 
est band with both ^ and ^ particles the system becomes 
a quantum spin Hall insulator with opposite chirality for 



Figure 3: (color online) Artificial magnetic fluxes 4> through 
the plaquettes of a square lattice (lattice constant indicated 
by the grid) resulting from superlattice modulation. Stripes 
or larger patches with strong, rectified magnetic fluxes can be 
achieved, (a) and (b): Single-component superlattices with 
(a) q x = (7v/d)(e x + e y /2), Vi = 4hu\ (b) q x = ±(7r/d)(e x + 
By/ 2) and V\ — 20hu. (c): Two components with q 1 / 2 = 
(ir/d) (e x =b e y ) , V1/2 = 12fiw. (d): Like (c), but with wave- 
lengths and amplitudes doubled. Always (p s = 0. 



the two species |36| . 

As an example, we consider unidirectional forcing 
F(t) = F(t)ep, with ep — cos((^ j p)e a; + sm((p F )e y 
and F(t) as depicted in Fig. [l]i (with T 2 = T/2 and 
thuo = AE/2, see supplemental material [30J for an an- 
alytical expression of the resulting phases). By varying 
the angle (pp and the forcing strength a = dFoTi/(27rh) 
(with lattice constant d) , we can access various topologi- 
cal quantum phase transitions, where at least one of the 
masses vanishes and changes sign (Fig. ^p). Thus, the 
lowest band can acquire a non-trivial Chern number. The 
inset shows how Dirac points can be moved and merged. 

A way to measure the topological band structure of 
the system is given by the method of Ref. [37J based on 
semi-classical wave-packet dynamics. It can be applied 
thanks to the adiabatic principle for Floquet systems [38J 
(see [39] for its application to the effective Hamiltonian). 

Superlattice modulation and flux rectification. In lat- 
tices with pairwise parallel bonds, such as square lat- 
tices, homogeneous driving vf(t) = —F(t) • as consid- 
ered in the previous paragraphs cannot create magnetic 
fluxes. Therefore, we propose to drive the system via 
an oscillating superlattice potential Vi(t) = f(t)Vo(ri) = 
f(t) J2 S if cos (^s ' r ~ where Vo(r) may be incom- 
mensurate with the host lattice. The driving function 
f(t) = f(t + T) breaks symmetries (r) and (s). To 
achieve a vanishing mean, (f(t))T — 0, in an experiment 
one can use 7r-shifted non- interfering standing waves such 
that f(t)V s cos(q s • r - <p 3 ) = V/(t)cos(q s • r - cp s ) + 
Vjf{t) cos(q s ■ r - tp 8 + tt), with V/, V/' > 0. In Fig. [3) we 
show - on the example of a square lattice with a shak- 
ing function as in Fig. [TJi (with T x /T = 0.8) - that, 
using different superlattice structures, various configura- 
tions of plaquette fluxes can be engineered [40]. Roughly, 
the larger the superlattice wavelengths the slower is the 
variation of the artificial flux. Therefore, superlattice 
modulation can generate not only strong magnetic fluxes 
through square plaquettes, but also large regions (stripes 
or patches) with rectified magnetic field where strong- 
field quantum Hall-type physics can be studied. Their 
inhomogeneity and finite extent provide a promising test 
ground for the investigation of robust edge modes. 
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Figure 4: (color online) Non-Abelian SU(2) gauge fields, (a) 
Two standing laser waves (with a phase shift of 71/ 2 and in- 
plane polarization as denoted in the figure) create a bipartite 
square lattice with alternating cr + and a~ polarized sites (A 
and B) |41| . rriF = ±1 particles feel an energy difference of 
±AE between A and B sites, (b) The resulting level scheme. 
A constant B-field realizes an additional on-site energy split- 
ting AE' (green arrow) such that \AE a ,b\ = | ± AE + AE'\ 
becomes sublattice dependent. The coupling Q of both spin 
states can be realized by magnetic or microwave fields, (c) 
Trace of the Wilson loop L in parameter space. Deviations 
from 2 imply non-Abelian physics [K y = 1.814; outside the 
white (black) regions, |trL| < 1.9 (< 1.99)]. 

Non-A belian SU(2) gauge fields. The periodic driving 
also permits the creation of arbitrary non-Abelian SU(2) 
gauge fields. Consider t and \, particles (say, mp = ±1) 
loaded into the spin-dependent square lattice depicted in 
Fig. [4^i, where the energy of t particles is lifted (lowered) 
by AE/2 on A (B) sites, and vice versa for ^ particles. 
These energy shifts are summarized by AEa z s z /2, if we 
introduce two sets of Pauli matrices s and a for spin 
(t* or 1) and sublattice (A or £?), respectively. More- 
over, uniform microwave and magnetic fields can be em- 
ployed to couple the f and | state with a matrix ele- 
ment ft and to produce an additional site-independent 
energy splitting AE' ', giving the site-independent term 
AE'sz/2 + Vts x . The absolute value of the total 
splitting AEi = AEa z — AE' is sublattice-dependent 
(Fig. [IJd). Including the NN tunneling J and a spin- 
independent sinusoidal drive vf(t) = —ri ■ Focos(ujt) 
as it can be induced by simply shaking the lattice back 
and forth, the Hamiltonian reads H = J^l^j + 

J2 i al[lAE i s z ^ns x ^vf(t)]d h with at = (a^aj^). The 

transformation bi = u\cbi, where ui are time- independent 
unitary 2x2-matrices, diagonal izes the Ha miltonian on 
site with eigenvalues h\ = | ^ AEf + 4Q 2 . This yields 

H = -E W) JbUtujbj + EMtiXiS z + vnWi. The 
sublattice dependence of Ui through AEi/ (20) achieves 
generically u\uj ^ 1. As in the derivation preceding 
Eq. |l]), the unitary transformation exp ( — i J2 i b\[Xits z — 
Ki sm(ujt)]bij with K{ = Vi ■ Fo/(huj) leads to a purely 
kinetic Hamiltonian H' = — J2(ij) Jb\Wij(t)bj. Here, 

W.. m _ e -i^sinM) ( c^^y dij e-<^i>\ 



and Cij and parametrize u\uj. For Huj ^> Jij, we 
can approximate H' by its time average H e ^ = {H')t — 
~^2(ij) JiffilMijbj , with the effective tunneling ma- 
trix elements JfJ* = J^/j det( (W^t)!, and the matrices 

Mij = {W i3 ) T /^\det({W i3 )T)V For J5f ^ 0, we require 
\ e B ± A^ e A = z/± Co> with integers ^± , and for unitarity of 
, we require v± both either odd or even. 

If the so called Wilson loop L, the product of the matri- 
ces M^ around a plaquette, yields not just a simple phase 
e^H describing an Abelian magnetic flux </>, the system 
is subjected to a genuine non-Abelian SU(2) gauge field. 
This is equivalent to requiring \tiL\ < 2, a sine qua non 
for the anomalous integer quantum Hall effect [42] and 
fractional quantum Hall states with non-Abelian anyonic 
excitations [43J. Without driving, |trL| = 2, but includ- 
ing it, |trL| < 2 can be fulfilled (see [3D]). 

Let us choose v + = 3 and v_ = 1, achieved by 

AE B = a/ 4 (M 2 + ae a and n = \/(M 2 ~ ae a/^. 
This leaves AEa/^, K x , and K y as free parameters 
(where K XjV is the amplitude of the forcing Kij in positive 
x, ^/-direction). In Fig. [4]:, we plot the trace of the Wil- 
son loop |trL| versus K x /K y and Ea/TujU for K y = 1.84 
(this value is not crucial but ensures large y-tunneling 
- see supplemental material [30], where also an analyti- 
cal expression for the Wilson loop is derived). There are 
broad regions where |trL| differs strongly from 2, prov- 
ing the presence of a strong artificial non-Abelian gauge 
field. Under typical conditions, the system shows Dirac 
cones, be it Abelian or non-Abelian. Similar analytic cal- 
culations reveal that L = t in a hexagonal lattice (and 
for even v± in the square lattice) . This limitation can be 
overcome by employing position-dependent coupling via 
Raman laser mixing, Q — »> Qi = Qe iq Vi with q the laser 
wave- vector difference (see [3D]). This way, the M^ as 
well as L can be tuned to be a generic (ix) SU(2) ma- 
trix both in square and hexagonal lattices. Alternatively, 
in a hexagonal lattice a non-trivial Wilson loop can be 
achieved with NNN tunneling. 

Conclusion. The creation of artificial Abelian and 
non-Abelian gauge fields by means of time-periodic forc- 
ing opens realistic perspectives for experimental stud- 
ies. This method offers great flexibility, because it does 
not involve the internal atomic structure. For fermions, 
where only different internal states interact with each 
other, this can be very advantageous for reaching the 
strongly correlated regime. 
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Supplemental Material 



I. ANALYTIC FORM OF THE EFFECTIVE 
TUNNELING AND THE PEIERLS PHASES 



In the examples of the main text, we use a driving 
potential F(t) = F(t)eF which is unidirectional, i.e., 
ep = cos((pF)e x -\-sm(cpF)e y , and has a paused- sine- wave 
amplitude as depicted in Fig. Id, 



F(t) 



_ J F sin(27rt/Ti), 
0, 







t mod T < 
tmodT < 



T 



(SI) 

Carrying out the time integrations as given in the in- 
troduction of the main text, this driving creates - for a 
vanishing energy difference between sites i and j, char- 
acterized by Uij — - the effective tunneling 



II. NON-ABELIAN SU(2) GAUGE FIELDS 

A. Independence of the observables on the choice 
of the phases of the local basis 

In the main text, it is shown how a non-Abelian gauge 
field can be induced via lattice shaking in a bipartite, 
spin-dependent lattice. Here, we comment on some tech- 
nical and background aspects that have been omitted 
there. We adopt the same notation as in the main text. 
First, we focus on the the example of site-independent 
magnetic mixing = O, Mi. The on-site Hamiltonian is 



Hi = EjS z 



Qs x 



(S6) 



Teff rp 

ij _ i l ( 

J a T 



where Ei = Ea, if the site i belongs to sublattice A, while 
E{ = Eb-, if it belongs to sublattice B. The transforma- 
tion Ui diagonalizing Hi can be written as 



and for ^ the AC-induced tunneling (ACT) 



(S2) 



(S7) 



eff 



T 



llXVi 



^^JvtMj) (S3) 



Here, is the BesselJ function of order /i, 7^ corre- 
spond to the freedom to choose the local phases, and 
we defined the dimensionless driving amplitude ol^ = 



where = atan^. The phases (pi and Ti correspond to 
the jf : they attest the freedom in choosing the phases of 
the two states which form the local basis. As these are 
arbitrary, physical observables - for instance the Wil- 
son Loop operator L - cannot depend on such phases. 
This is immediate in the absence of periodic driving. In- 
deed, in that case the Wilson loop is the identity what- 
ever choice of Ui is performed, as it corresponds to prod- 
ucts of terms uiu\ = 1, for any site i in the loop. 
In presence of periodic driving, the cancellation of the 
phases is slightly more involved. Following the definition 



r ij -e F FoT 1 J(27rh). For the case T x = T/2 (and choosing of the effective hopping matrices M i3 in the main text 

Mij = (Wij)T / y/\deb((Wij)T)\, we notice that that their 



the local phases 7^ = 0), as considered in the section 
Topological and quantum spin Hall insulator, the Peierls 
phases Oij are thus for nearest neighbors (NN) (where 



matrix elements have the same phases (up to multiples of 



i) 



tt) as the elements 



while their moduli 



fiTij) = 2 (n- a ij) i 



(S4) 



L l U 3 ~ (-d% C *j)> 

are independent of such phases. This implies that for two 
different choices of the local phases at the site fc, ip f k , r' k , 
and (fk, Tk, respectively, the effective hopping matrices 
relate as 



and for next-NN (NNN) (where = 0) 



M' 

13 



i(A(fjS z -\-Arj 



(S8) 



tan 



OLi 



%^tan/(r,) = l±||^t a n ^ 



(S5) 



The flux threaded through a triangular pla- 
quette as sketched in Fig. la-b is then <p = 
f(e x ) + f(-e x /2 + V3/2e y ) + f(-e x /2 - V3/2e y ) 
(and similarly for Fig. lc, where the plaquette is 
spanned by NNN tunneling). It can be non-zero since 
f(rij) is a non-linear function of r^-. 



where Acpk = <p f k — (fk, an d ^ T k = T k — r^. That is, 
the choice of the phases commutes with the time-average 
procedure. It follows that they cancel out when the hop- 
ping matrices are multiplied in the Wilson Loop L as 
in the time-independent case. Hence, L is independent 
of the choice of local phases, as it should be. 

We conclude by remarking that the same happens for 
more involved choices of the optical lattice, as the actual 
form of the local Hamiltonian does not play any role (cf. 
section II C). 
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B. Analytic calculation of the Wilson Loop 



Here, we show how to derive analytically the Wil- 
son Loop L computed along the fundamental plaquette 
for a forced lattice described by (S6). Using the re- 



sult explained in the previous section, we may choose 
tpi = Ti = 0, Vi. Using this in (S7), this implies that 



u\uj = e2 



(Aj-A^Sy 



(S9) 



We may call S the forcing-induced non-linear map that 
relates 



• j to the effective hopping matrix , — 
Because we employed the unitary transforma- 
tion exp ( — i^Zibll^itSz — Ki sm(ujt)]bi) to arrive at 
Eq. (2) of the main text, this map is a function of the 
energy differences between A and B sites, A^b ± KeA = 
v±uj. To ensure finite ACT, the v± have to be inte- 
gers. Explicitly, as defined in the main text, = 
(Wij)T/V\ det (( W ij)T)l which yields, following Eq. (2), 



(Wij)T 



(cijJ^iKij) d i3 J u+ 

u 



Jv 



)■ 



(S10) 



Here, the amplitude of the forcing K{, 



F /(huj) depends only on the direction 5/ of link ij 
(81=1,2 = x,y for the square lattice, and ^=1,2,3 = 
x,—\{x — \^3y) : — |(x + V3y) for the hexagonal one). 



To ensure unitarity MijMj, 



1, we therefore have 



to impose the additional condition that the v± are ei- 
ther both even or both odd (use V-(ji) = —V-(ij) and 

First, we characterize the action of S on Eq. ( |S9| ) for 
even forcing, i.e., when both v+ and V- are even num- 
bers, v+ — 2n, V- — 2n' . By the observation that S 
maps by construction i) unitary matrices to unitary ma- 
trices, and ii) real matrices to real matrices, it follows 



that S[e i(fi *y] = e i(fi Sy , where <p' = atan[ 



tan (p] . 



The main consequence of this form of S[e ly:,Sy ] is that 
hopping matrices M^- = Mi commute also for different 
links /, and this implies that the Wilson loop for the 
square, L = MiM2MlM^ and for the hexagonal lattice 
L M 1 MlM2M\M ?) M\, are both equal to the identity 
(by convention, we travel the loop on fundamental cells 
anticlockwise, starting from the bottom-left corner). 

The situation is more interesting for odd forcing, 
i.e., when both v+ and v_ are odd numbers, v+ = 
2n + 1, V- = 2n' + 1. The result of the time- 
averaging is in this case S[e %ipSy ] = s z e Z(p Sy , where ip' = 
atan [ '^j^+iIk-) ta* 1 ^]- Due to the presence of the s z , for 
the square lattice, starting by convention from sublattice 
A, we get 



m x m 2 m\m\ 



s z e liPxSy s z e^ y ° y e 



s z e 



= e i( <- 



where ip r = atan 



tan [| (atan ^ atan^-)] 



H = x,y. As shown numerically in the main text, the 




Figure SI: Bloch sphere representation of the Wilson loop L 
computed for odd forcing in a hexagonal lattice with \q\ = 1 
(in lattice step units). Eb is taken to be Eb — — yj 4:E A + 3Q 
such that the eigenvalues enjoy the relation Xb = 2Xa = 
\/E\ + Vi — 2uj, which implies Xb + A a — 3cj, and Xb — Xa — 
uj. Here, uj is the frequency of the periodic forcing. The 
parameters q — (cos# q , sin# q ), the intensity of the shaking 
K = \K\(cos0k, sin Or), and Ea, are chosen randomly in the 
ranges, q G [0,2tt[, |K| G]0,5[, 6>k G [0,2tt[, and G]5, 15[, 
respectively. The resulting L's cover the entire Bloch sphere. 



difference cp' — <p' x can take any value between zero 
and 27r, allowing for a non-trivial Wilson Loop, 
tr|L| = 2 |cos(<^ — (p x )\ 7^ 2. In the example, we have 



chosen K,, 



1.84, because this lies close to the first 



maximum of the BesselJ function J\(K y \ ensuring large 
tunneling J?? = Jyf \ deUjWij)r)\ in ^/-direction. 

On the other hand, for the hexagonal lattice we get 
always a trivial result as the s^'s cancel out and L = 1. 



C. Raman- induced mixing 

Considering a site-independent mixing of |- and |- 
particles as in (S6), the form of Mij is limited to e lip Sy 
or s z e tipSy . This can be circumvented by considering 

In this case, the local 



site-dependent Raman mixing. 
Hamiltonian takes the form 



Hi = EiS z + ft [cos(q • n)s x - sin(q • ri)s y ] , (Sll) 

where q is the difference of the wave vector of the two 
Raman lasers, E { = E A , Mi G A, and Ei = E B , Vi G B. 

The main difference in this case is that, while still de- 
pending only on the link direction — rj, the M^ for 
different link directions are generally not commuting. In- 
deed, in this case the local transformation U{ may be cho- 
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sen as m — e" 



z e2 c 



j and all its neighbors z, the products i^i^- for differ- 
ent link directions are not commuting, and the map S 
acts highly non-trivial on them. Numerical studies with 



'. Hence, considering a site random selected parameters (Fig. SI) show that the cor- 



responding Wilson Loop L is in general non-trivial and 
dense in the Bloch sphere, i.e., we can simulate any de- 
signed L G SU(2), not just matrices of the form e l6Sy . 



